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Black holes generate collimated, relativistic jets which have been observed in 
gamma-ray bursts (GRBs), microquasars, and at the center of some galaxies 
(active galactic nuclei; AGN). How jet physics scales from stellar black holes in 
GRBs to the supermassive ones in AGNs is still unknown. Here we show that 
jets produced by AGNs and GRBs exhibit the same correlation between the 
kinetic power carried by accelerated particles and the 7-ray luminosity, with 
AGNs and GRBs lying at the low and high-luminosity ends, respectively, of the 
correlation. This result implies that the efficiency of energy dissipation in jets 
produced in black hole systems is similar over 10 orders of magnitude in jet 
power, establishing a physical analogy between AGN and GRBs. 

Relativistic jets are ubiquitous in the cosmos and have been observed in a diverse range of 
black hole systems spanning from stellar mass (~ 10 M Q ) to supermassive scales (~ 10 5 — 
10 lo M e ), in particular in the bright flashes of gamma-rays known as GRBs ( 1 , 2 ), the miniature 
versions of quasars lurking in our galaxy known as “microquasars” (5) and AGNs ( 4 , 5 ). Despite 
decades of observations at almost all wavelengths and considerable theoretical efforts, there are 
still many aspects of black hole jets which remain mysterious: the mechanism(s) responsible for 
their formation and the nature of their energetics as well as their high-energy radiation ( 6 , 7 ). Jets 
and outflows from supermassive black holes have important feedback effects on scales ranging 
from their host galaxies to groups and clusters of galaxies (8). Hence, a better understanding 
of the physics of jets is required in order to have a more complete picture of the formation 
and evolution of large-scale structures in the universe and the coevolution of black holes and 
galaxies (9). 

One outstanding question is how the jet physics scale with mass from stellar to supermassive 
black holes. Interestingly, there is evidence that jets behave in similar ways in microquasars and 
radio-loud AGN ( 10 - 12 ). However, a clear connection between AGN and GRBs has not been 
established yet, although recent work provides encouraging results ( 13 , 14 ). 

As a first step in understanding how the properties of jets vary across the mass scale, we 
focus on the energetics of jets produced in AGNs and GRBs. Therefore, we searched the litera- 
ture for published and archival observations that allow us to estimate the jet radiative output and 
the kinetic power for a sample of black hole systems in which the jet is closely aligned with our 
line of sight and characterized by a broad range of masses. For this reason, our sample consists 
of blazars - AGNs with their jets oriented toward Earth ( 15 ) - and GRBs, the spectral energy 
distributions of which are completely dominated by the jet due to beaming effects. 

We used as a proxy of the jet bolometric luminosity the observed 7-ray luminosity L lso 
which is isotropically equivalent. In order to estimate the kinetic power Pj et , we use extended 
radio luminosities for the blazars whereas for the GRBs we relied on the afterglow measure- 
ments in radio or X-rays. Therefore, the availability of these observables restricted our sample 
to 234 blazars (106 BL Lacs and 128 flat-spectrum radio quasars - FSRQs; see Table SI) and 
54 GRBs (49 long and 5 short GRBs, all with known redshifts z\ see Table S2). For blazars, 
L lso was estimated from the 7-ray energy flux and the spectral index measured with Fermi Large 
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Area Telescope (LAT) (76); P- ]et was estimated using an empirical correlation which relates the 
Very Large Array (VLA) extended radio emission and the jet kinetic power (77, 18). For GRBs, 
L lso = E lso ( 1 + z)/tgo where t< M is the burst duration and P ISO is the isotropically equivalent 
energy radiated during the prompt emission phase and measured with different telescopes (21 
observed with either BeppoSAX, BATSE, HETE, HETE-2 or Integral, 24 with Swift Burst Alert 
Telescope - BAT - and 10 with Fermi). Pj et was computed as Pjet = fbE%°( 1 + z)/t 90 where 
7% so is the kinetic energy estimated from the radio (VLA) or X-ray (Chandra) luminosity during 
the afterglow phase using the standard afterglow model (79), f b = 1 — cos 0 is the “beaming 
factor” and 9 is the radiation cone half-opening angle which is the same as the jet opening angle 
estimated from the jet break in the GRB afterglow lightcurve (20). 

We first compared the relative trends of L lso and ^jet for the blazar and GRB population 
separately (Figure 1). The Pearson correlation coefficients of 0.85 and 0.8 obtained for blazars 
and GRBs respectively, indicate a strong correlation within each group of sources. However, 
the L is °-Pj et trend is different for GRBs and blazars as shown by the fits to the data (Fig. 1). 

We computed the intrinsic luminosity, L, for GRBs and blazars by correcting L lso for the 
opening angle or beaming factor, f b , such that L = fbL lso . For GRBs, the beaming factor 
is computed from the jet opening angle Oj as 1 — cos Oj (22); for blazars, f b is estimated as 
1 — cos 1 /r where T is the bulk Lorentz factor of the flow, since AGNs obey 7, < 1/T (23, 24). 
While an estimate of Oj is available for each GRB in the sample, F is only available for a subset 
of 41 blazars. Figure 2 shows an anti correlation between L lso and f b for both GRBs and blazars 
with compatible indices when fit with a power law. Because 9 is not available for the whole 
blazar sample, we used the power-law fit of L lso vs f b as an estimator for f b . 

As with L lso and P jet , L and Pj et are strongly correlated within the GRB and AGN samples 
(Fig. 3). However, they follow the same trend within the narrow uncertainties and the whole 
GRB and blazar sample can be fit adequately with a power law over 10 orders of magnitude 
in luminosity. Therefore, the relativistic jets in GRBs and blazars are consistent with obeying 
the relation P je t « 4.6 x 10 47 (L/IO 47 ) 0 ' 98 erg s -1 , within the measurement uncertainties. In 
other words, once “black hole engines” produce relativistic jets, they seem to do so maintaining 
the same coupling between the total power carried by the jet and power radiated away. This 
universal scaling for the energetics of jets is maintained across the mass scale regardless of the 
different environments and accretion flow conditions around the compact object. 

Figure 4 indicates that most of the jets in our sample dissipate at least 3% of the power 
carried by the jet as radiation and overall they can radiate as much 15%. This range of ef- 
ficiencies is considerably higher than previous estimates for AGNs based on radio to X-rays 
luminosities (25, 26) but they are in agreement with results obtained from blazar broadband 
spectral models (27, 28) as well as GRB afterglow studies (29-31). Efficient heating of elec- 
trons seems to be a universal property of relativistic magnetized shocks according to numerical 
simulations (32) which demonstrate that electrons retain > 15% of the pre-shock energy. If 
most of the post-shock energy is radiated away, these theoretical results could pave the way to 
an understanding of the high dissipation efficiencies that we find. 

Our results suggest that there is a single fundamental mechanism to produce relativistic jets 
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in the Universe. The analogy known to exist between microquasars and AGNs ( 3 , 10 , 11 ) can 
be extended to the gamma-ray bursts with the fundamental difference that whereas AGNs and 
microquasars undergo recurrent activity, GRBs experience only one episode of hyperaccretion. 
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Figure 1: The relation between the jet kinetic power and the isotropically-equivalent 7-ray 
luminosity for AGNs and GRBs. Error bars, lcr. We fitted the two populations separately 
using a symmetric least-squares method (orthogonal BCES with bootstrapping; 21). The blazar 
and GRB best-fit models corresponds to the solid and dashed lines, respectively (log P jet = 
A log L lso + B ). The best-fit parameters obtained for the blazars are A = 0.51 ±0.02 and 
B = 21.2 ± 1.1; for the GRBs, A = 0.74 ± 0.08 and B = 11.8 ± 4.1. The scatter about the 
best-fit is 0.5 dex and 0.8 dex for the blazars and GRBs, respectively. The 2<r confidence band of 
the fits is shown as the gray shaded regions (barely visible for blazars). The two correlations do 
not agree at > 5cr level. We also include for illustration XRF 020903 and GRB 090423 (yellow 
circles) as well as the two recent tidal disruption flares (TDFs) detected with Swift which are 
presumably due to the onset of relativistic jets from the tidal disruption of stars by supermassive 
black holes (36). We do not consider these sources in the statistics since we only have limits on 
their luminosities 
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Figure 2: The relation between the apparent 7 -ray luminosity and the beaming factor for blazars 
(left panel) and GRBs (right panel). We find r = —0.53 and —0.56 for blazars and GRBs, 
respectively, indicating anti correlations significant at the 3.6<r and 4.4<r levels respectively. The 
solid lines correspond to the best-fit linear models obtained with the symmetric least-squares fit 
and are given by /& « 5 x 10 -4 (L4 S 9 )-°- 39±0 - 15 and ~ (J.(J3(L)g’) _0 ' 24±0 0fj for blazars and GRBs, 
respectively. The gray shaded region corresponds to the lcr confidence band and the blue and 
yellow regions are the la prediction bands, which quantify the scatter about the best-fits. 
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Figure 3: The relation between the collimation-corrected 7-ray luminosity L = ft,L lso and the 
kinetic power for AGNs and GRBs. The shaded regions display the 2 a confidence band of 
the fits. The blazar and GRB best-fit models (dashed and dotted lines, respectively) follow 
correlations which are consistent, within the uncertainties, with the best-fit model obtained 
from the joint data set (solid line). In other words, using L instead of L lso leads to correlations 
for AGNs and GRBs which are consistent with each other (compare to Fig. 1). The best-fit 
parameters obtained from the combined data set are a = 0.98 ± 0.02 and f3 = 1.6 ± 0.9 
where log Pj et = a log L + j3. The scatter about the best-fit is 0.64 dex. The yellow data 
points correspond to XRF 020903 and GRB 090423, which we do not take into account in the 
statistics. 
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Figure 4: Distribution of the la lower limits on the jet radiative efficiency (the fraction of the 
total jet power which is converted to 7 -rays) e ra d = L/(L + Pj et ) for AGNs and GRBs. The 
vertical solid lines indicate the median values of the lower limits and the dashed lines represent 
the median values of e ra d for each sample. Most of the sources are characterized by e ra d > 3%. 
The median efficiencies correspond to about 15%, keeping in mind that these estimates are 
affected by ~ 0.5 — 0.7 dex uncertainties on average. 
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Supplementary Materials 

Data 

Table SI lists the properties of the 234 blazars in our sample (106 BL Lac objects and 128 
FSRQs). Our sample includes the well-studied blazars OJ 287, 3C 454.3, 3C 279 and 3C 273. 
The redshifts were estimated as described in (77). 

Table S2 lists the corresponding properties of the 54 GRBs (49 long and 5 short). We 
include in our data the sub-energetic bursts GRB 031203 (37) and GRB 980425 (38). The latter 
one is associated with the nearby (distance ~ 40 Mpc) supernova 1998bw (39, 40). We also 
include in our sample the naked-eye GRB 0803 19B (42). We include in Table S2 the X-ray 
flash (XRF) 020903 (41) and the most distant cosmic explosion ever detected, GRB 090423 at 
z & 8.2 (43, 44) but we do not take these GRBs into account in the statistics, since we have only 
limits on their collimation-corrected energetics. 


Calculation of 7-ray luminosity and jet power 

Blazars 


7 -ray luminosity: 

In order to calculate the total Fermi 7-ray luminosity, we follow a procedure similar to that 
of (46) (their equation 1). The procedure to calculate the /[‘-corrected band luminosity depends 
on the type of model used in the 2FGL analysis (16,47). In the case of a power law energy 
model, the total 100 MeV to 100 GeV /c-corrected luminosity is calculated from the energy flux 
(5' 7 ) given in the catalog: 


L iso 


4ircP L 


ff 7 

(1 + zy~ a ’ 


(i) 


where d L is the luminosity distance in cm 2 and o 7 is the (energy) spectral slope over the whole 
band. 

In the 2FGL catalog, some blazars are now modeled with the ’Log Parabolic’ form. The 
^-corrected energy flux in this case must be calculated numerically. The integral form is 


s; = x 
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Eq(1 + z) 
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B 0 (l+z) 


(1 + z)~ 2 EdE, 
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where E\ = 0.1 GeV and E- 2 = 100 GeV and we have used the fit given in the 2FGL catalog for 
each source, with values a (column name ‘spectraLindex’), f3 (‘beta’), and E 0 (‘pivot.energy’), 
and K (‘flux_density’). The constant % = 1.6 erg MeV GeV -2 gives the energy flux in final units 
of erg cm -2 s -1 . 

The band luminosity for the Log-Parabolic model case can then be simply calculated from 


L iso = 47 Tdl s' r 


(3) 
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We calculate the uncertainty in L lso propagating the error associated with a 7 and S 1 quoted 
in the 2FGL. The average uncertainty in L lso corresponds to 0.05 dex. 

In order to estimate the uncertainty affecting the collimation-corrected luminosities L we 
first evaluate the error in the beaming angle 9 or correspondingly T. The uncertainty in T is 
dominated by the uncertainty in the variability Doppler factor whereas the uncertainty in the 
apparent speed does not contribute significantly to the error budget of T. The uncertainty in 
the variability Doppler factor is ~ 27% (1 s.d.; 48). Therefore, for the blazars with direct 
estimates of T available, the relative uncertainty in T is 0.3 (24, 48) which translates to an 
average uncertainty of 0.26 dex in L for these blazars. For the blazars without direct estimates 
of T, we estimate the uncertainty in L using the prediction band of the L lso — f b relation shown 
in Fig. 2. The plotted prediction band corresponds to a relative uncertainty of 0.69 in 9. The 
resulting average uncertainty affecting L for the blazars without direct estimates of Y is then 0.6 
dex. 

Kinetic power: 

Following (17), we estimate the jet kinetic power by using the correlation between the ex- 
tended radio emission and the jet power (18, 49). Cavagnolo et al. searched for X-ray cavities in 
different systems including giant elliptical galaxies and cD galaxies and estimated the jet power 
required to inflate these cavities or bubbles, obtaining the tight correlation 


between the “cavity” power and the radio luminosity. Hence, assuming Pj et = -P cav we can 
estimate the jet kinetic power for the blazars which have extended radio emission observed with 
the VLA (17). 

The uncertainty in Pj et is dominated by the scatter in the correlation of (18) and corresponds 
to 0.7 dex. 

GRBs 

7 -ray luminosity: 

The current scenario for GRBs (50, 51) posits that initially most of the energy produced by 
the GRB is in kinetic form produced during the short “active” state of the stellar-mass central 
engine. A certain fraction of initial energy is converted after a few seconds mostly to 7 -rays 
observed during the prompt emission, by means of internal shocks in the jet (52). The ultrarel- 
ativistic jet produced in the explosion later on collides with the circumburst medium producing 
the afterglow. Two crucial quantities which we use in this work are the radiative and kinetic 
energies released by the GRBs during their short period of activity. 

The isotropically equivalent energy radiated in 7 -rays E™° is directly available from mea- 
surements. It was measured for the GRBs using a variety of different telescopes including 
pre-Swift telescopes (BeppoSAX, BATSE, HETE, HETE-2 and Integral) as well as Swift and 
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erg s 


-1 
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Fermi (see Table S2). We calculate the isotropically-equivalent 7 -ray luminosity as 


L iso 


( 1 +*) 

^90 


E*° 


( 5 ) 


where t 90 is the duration containing 90% of the fluence in the observer frame. 

The energy range in which the fluence is measured is typically ~ 10 keV - 10 MeV. Most 
of the radiative energy released in the GRB jet during the prompt emission is contained in this 
energy range according to the latest GRB SEDs observed (55). 

We adopt an uncertainty of 0.2 dex on the values of Ef° which corresponds to the typical un- 
certainty affecting the GRBs in the sample studied by (54). Therefore, the resulting uncertainty 
in the value of L lso corresponds to 0.2 dex. 

The collimation-corrected 7 -ray luminosity is computed as 


L 


(1 + *) 

1 90 


f b Ef° 


( 6 ) 


where fb is the beaming factor (J), = 1 — cos f )- } ) and 9 j is the jet half-opening angle. This relies 
on the afterglow lightcurve displaying a jet break which is used to estimate 9 (22). (9j) ~ 8 ° 
for the sample and (fb) ~ 9 x 10 -3 . We adopt an uncertainty of 0.1 dex on the values of #jet; 
which corresponds to the typical uncertainty in the values of 9 j for the sample studied by (54). 

We calculated the uncertainty in L using error propagation from the uncertainties in Ef° 
and 9, obtaining that the uncertainty in L is ~ 0.3 dex. 

Kinetic energy: 

The jet kinetic energy is estimated from the radio or X-ray afterglow lightcurve using the 
fireball model (19). For most of the GRBs, X-ray data were used to determine this energy 
(31, 54). In a few cases, radio data were used (35, 37, 41). 

The standard fireball afterglow model depends on five model parameters: the explosion 
kinetic energy Ef°, the density of the circumburst environment n (with which the jet collides), 
the spectral index of the electron energy distribution p and the fractions e e and eg of shock 
thermal energy carried by electrons and magnetic field, respectively. The typical values adopted 
for these parameters are p = 2.2, n = 1 cm -3 , e e = 0.1 and e B = 0.01 (e.g., 19, 31). The 
afterglow model relates the specific flux at a certain frequency and at a specific time after the 
burst (typically 10 hours in the observed frame) with the kinetic energy. 

The measurement of E\ f° can be impacted by the afterglow plateau which possibly corre- 
sponds to a late activity of the central engine. Indeed, (30) demonstrated the impact of choosing 
two different times for the measurement of i% so : Ticc (deceleration time) or f 6 (injection break 
time). Adopting t dec would lead to an underestimation of Ef° because this choice of time does 
not include the plateau. 

The kinetic energy estimates we used in our analysis were computed at either ~ 10 h or 
~ 24 h (cf. Table S2). We verified that these times are usually beyond the “break time” of the 
plateau reported in (45). Therefore, the measurements of Ef° used in this paper correspond to 
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conservative estimates. Moreover, the choice of the time at either ~ 10 h and ~ 24 h affects 
very little the measurement of ° (31). 

The typical uncertainty in Ef° due to observational errors is ~ 0.3 dex (30, 54). The value 
of Ef° is also sensitive to the values of parameters which regulate the microphysics of the 
fireball afterglow model and are poorly constrained. The systematic error affecting E™° due to 
the uncertainties in the parameters e e and e B of the fireball afterglow model can be as high as 
0.45 dex (29, 30). Therefore, we combined the uncertainty resulting from the observational and 
systematic sources of errors in quadrature and conservatively adopt an uncertainty of 0.5 dex 
for Ef°. 

Kinetic power: 

The jet kinetic power is computed as 

Ph = (7) 

t-90 

L™° and Pj et should be thought as the average luminosities over the duration t< M of the prompt 
emission phase, i.e. the average luminosities over the timescale during which the central engine 
is producing the jet. We calculated the uncertainty in ^jet using error propagation from the 
uncertainties in PJf’ and 0, obtaining that the uncertainty affecting P jet is ~ 0.54 dex. 

Linear regression method 

Here we present more details about the linear regression method that we use in the paper. 

We fitted the datasets using the BCES (bivariate correlated errors and intrinsic scatter) re- 
gression method (21) which takes into account measurement errors in both the “X” and “Y” 
coordinates and the intrinsic scatter in the data. This method has been widely used in fitting 
datasets in the astronomical community (70, 71). 

It is not clear in the data sets analyzed in this work which quantities should be treated as 
the dependent variables and which ones should be treated as independent from a physical point 
of view. There is no a priori reason to expect the luminosity to be the independent variable 
as opposed to the kinetic power (or beaming factor). For this reason, we treat the variables 
symmetrically and adopt the BCES orthogonal regression method, which minimizes the squared 
orthogonal distances. Uncertainties on the parameters derived from the fits are estimated after 
carrying out 100000 bootstrap resamples of the data. 

Partial correlation analysis 

When studying correlations between luminosities one should be careful to take into account 
their common dependence on the distance (11). We performed a partial correlation analysis of 
the common dependence of Ei)° and Pj et on the distance using the partial Kendall’s r correlation 
test (72). 
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Objects 

N 

T 

a 

-Pnull 

Signif. rejection null 

Blazars 

234 

0.3 

0.04 

5 x 10“ 15 

7.8 

GRBs 

54 

0.4 

0.08 

1.6 x 10“ 7 

5.2 


Table S3 Results of partial correlation analysis with X = Pj et , Y = L lso and Z = log d p 
Notes. Column (1): subsample. Column (2): Number of sources. Column (3)-(6): results of 
partial correlation analysis; r is the partial Kendall’s correlation coefficient; a is the square root 
of the calculated variance; P nu n is the probability for accepting the null hypothesis that there 
is no correlation between X and Y ; Column (6) gives the associated significance in standard 
deviations with which the null hypothesis is rejected. 


Objects 
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T 

a 

■Pnull 

Signif. rejection null 

All objects 

288 

0.57 

0.06 

<C 10- 10 

9.4 

Blazars 

234 

0.24 

0.03 

2 x 10- 12 

7 

GRBs 

54 

0.47 

0.08 

10“ 9 

6.1 


Table S4 Same as Table S3 with X = Pj et , Y = L and Z = log d L . 


We applied this test to our data considering X = P je t , Y = L lso and Z = log d L . Table S3 
lists the results of the partial correlation analysis. This test demonstrates that the p - value of the 
null hypothesis (i.e. no correlation between X and Y) is 1.2 x 10 -7 when considering the GRB 
subsample and < 10 -10 when considering the blazar subsample or the combined blazar-GRB 
sample. Therefore, the L lso — P jet correlation is strong and not a distance-driven artifact. 

We then applied the partial correlation test to the data considering X = Pj et , Y = L and 
Z = log d L , i.e. we replaced the isotropically-equivalent luminosity in Y with the collimation- 
corrected one. Table S4 lists the results of the analysis. This test demonstrates that the p- value 
of the null hypothesis is <C 10 _1 °, ~ 10 -12 and 10 -9 when considering the combined blazar- 
GRB sample, the blazars and the GRB, respectively. Hence, the L — P jet correlation remains 
very strong after correcting for beaming. 

Blazar luminosity estimates: Impact of the synchrotron peak 

We discuss in this section the impact of the lower energy synchrotron peak - observed in the 
spectral energy distribution of blazars (74, 75) - in estimating the jet radiative luminosity. 

A subset of 13 1 blazars - roughly half of the original AGN sample - have adequate sampling 
of their multiwavelength spectral energy distributions which allow us to quantify the impact of 
the synchrotron peak in the estimate of the jet radiative luminosity. We estimated the bolometric 
luminosity as Lj^j = L lso + Lf° n where Lf° n is the isotropically-equivalent luminosity of the 
synchrotron peak. We obtained that the values of Ljp for the quasars are not much different 
from L lso (on average by a factor of ~ 1.7), whereas the values of Ljp for the BL Lacs are 
somewhat different (on average by a factor of ~ 2.8). 
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We computed the intrinsic bolometric luminosity as L ho \ = correcting L™° l for the 

opening angle or beaming factor. As previously discussed, estimates of fb for blazars rely on 
measurements of T which are not available for all the blazars in our sample. For this reason, 
we use the anti correlation between Lj^ and fb as an estimator of fb for the blazars without 
measurements of T. 

We show in Figure SI the resulting relation between L ho \ and Pj Ct compared to the L — P- ]Ct 
fit derived before (cf. Fig. 3). Figure SI illustrates that the L bol — P jet and L — P jet best-fits are 
characterized by very similar slopes. However, the fit based on Lboi is characterized by slightly 
higher radiative efficiencies: for a given jet power, the jet luminosity is higher on average by a 
factor of « 2 compared to the fit based on L. Therefore, the results based on L bo i strengthen 
our conclusion that AGN jets have high radiative efficiencies while still being in qualitative 
agreement with the GRB result. 
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Figure SI The relation between the collimation-corrected bolometric luminosity L bo i and the 
kinetic power for the 131 blazars with measurements of synchrotron peak luminosity. The 
solid line in this figure corresponds to the fit based on L ho \ whereas the dashed line is the fit 
based on L. The shaded regions display the 2 a confidence band of the fits. Error bars, la. 
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Table SI. Properties of the blazar sample. 
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J2347.9-1629 PKS 2345-16 FSRQ 0.576 46.5 0.05 44.14 0.57 45.21 -2.36 

J01 16.0-1 134 PKS 0113-118 FSRQ 0.672 46.71 0.03 44.28 0.57 45.46 -2.43 

J2229.7-0832 PKS 2227-08 FSRQ 1.56 48.26 0.03 45.96 0.26 45.17 -2.3* 
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J1012.1+0631 BZB J1 01 2+0630 BLL 0.727 46.3 0.09 44.01 0.57 45.32 -2.29 

J1023. 6+3947 B2 1020+40 FSRQ 1.254 47.02 0.08 44.48 0.57 45.78 -2.54 

J1051. 3+3938 BZB J1051+3943 BLL 0.49 45.57 0.17 43.53 0.59 43.93 -2.04 

J1040.7+0614 4C +06.41 FSRQ 1.27 47.38 0.05 44.72 0.57 45.06 -2.66 


J 1058.4+0133 4C +01.28 BLL 0.888 47.42 0.02 45.23 0.26 45.61 -2.19* 

J1 104.4+3812 MKN421 BLL 0.03 44.88 0.01 43.08 0.56 43.28 -1.8 

J1 121.0+4211 EXO 1118.0+4228 BLL 0.124 44.82 0.07 43.04 0.57 42.61 -1.78 

J1 132.9+0033 PKS 1130+008 BLL 1.223 47.32 0.05 44.68 0.57 45.46 -2.64 

J1130. 3-1448 PKS 1127-145 FSRQ 1.187 47.68 0.03 44.92 0.57 45.62 -2.76 
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J1337.7-1257 PKS 1335-127 FSRQ 0.539 46.53 0.04 44.16 0.57 45.15 -2.37 

J1354.7-1047 PKS 1352-104 FSRQ 0.332 45.93 0.04 43.77 0.57 44.84 -2.16 

J 1405. 1+0405 PKS 1402+044 FSRQ 3.215 47.68 0.16 44.92 0.59 46.29 -2.76 

J1419.4+3820 B3 1417+385 FSRQ 1.82 47.58 0.09 44.85 0.57 45.07 -2.73 


J 1418. 1+2539 BZB J14 17+2543 BLL 0.237 44.92 0.12 43.1 0.57 44.2 -1.81 

J 1427. 0+2347 PG 1424+240 BLL 0.16 45.98 0.02 43.8 0.56 44.22 -2.18 

J1440. 3-1540 PKS 1437-153 BLL 0.244 45.1 0.1 43.22 0.57 44.59 -1.88 

J 15 13.6-3233 CRATES J15 13-3234 FSRQ 0.244 45.6 0.04 43.55 0.57 44.79 -2.05 

J1517.7-2421 APLib BLL 0.049 44.5 0.03 42.83 0.56 43.57 -1.67 
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J1 746.0+23 16 BZQ J1 745+2252 FSRQ 1.884 47.73 0.11 44.95 0.58 45.53 -2.78 

J1 8 1 3.5+3 143 B2 1811+31 BLL 0.117 44.85 0.04 43.06 0.56 43.62 -1.79 

J1 61 3.4+3409 OS 319 FSRQ 1.401 46.98 0.1 44.88 0.28 45.3 -2.11* 

J2004.5+7754 S5 2007+777 BLL 0.342 45.71 0.05 43.62 0.57 44.33 -2.09 


J1640.7+3945 NRA0 512 FSRQ 1.66 48.12 0.04 45.2 0.57 45.51 -2.91 

J1 635.2+38 10 4C +38.41 FSRQ 1.814 48.72 0.02 45.44 0.26 45.62 -3.29* 

J 1642.9+3949 3C 345 FSRQ 0.593 46.95 0.04 43.75 0.26 45.74 -3.2* 

J 1653. 9+3945 MKN501 BLL 0.034 44.46 0.02 42.8 0.56 43.2 -1.66 

J1749. 1+4323 B3 1747+433 BLL 0.229 45.38 0.05 43.41 0.57 44.03 -1.97 
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J1727. 1+4531 S4 1726+45 FSRQ 0.714 46.89 0.03 44.4 0.56 45.09 -2.49 

J1728.2+5015 IZW187 BLL 0.055 43.83 0.08 42.39 0.57 43.28 -1.44 

J1700.2+6831 BZQ J1700+6830 FSRQ 0.301 46.2 0.02 43.95 0.56 44.37 -2.25 

J0721.9+7120 S5 0716+714 BLL 0.3 46.78 0.01 44.46 0.26 44.76 -2.32* 


J07 10.5+5908 EXO 0706.1+5913 BLL 0.125 44.7 0.09 42.96 0.57 43.8 -1.74 

J0654.5+5043 CRATES J0654+5042 FSRQ 0.135 45.17 0.03 43.27 0.56 43.62 -1.9 

J0841. 6+7052 4C +71.07 FSRQ 2.172 48.36 0.05 45.16 0.26 45.97 -3.2* 

J0654.2+4514 S4 0650+45 FSRQ 0.933 47.34 0.03 44.7 0.56 45.1 -2.65 

J0747.7+4501 B 3 0745+453 FSRQ 0.192 44.78 0.1 43.02 0.57 43.87 -1.77 
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J1 15 1 .5+5857 BZB J1 151+5859 BLL 0.127 44.59 0.08 42.89 0.57 43.9 -1.7 

J1058. 6+5628 BZB J1058+5628 BLL 0.144 45.47 0.03 43.47 0.56 43.68 -2.0 

J1037. 6+5712 BZB J1037+5711 BLL 0.106 44.97 0.04 43.14 0.56 43.34 -1.83 

J0929.5+5009 BZB J0929+5013 BLL 0.37 45.57 0.08 43.53 0.57 44.1 -2.04 


J08 18.2+4223 S4 0814+425 BLL 0.245 46.17 0.02 45.01 0.26 44.41 -1.17* 

J1253. 1+5302 S4 1250+53 BLL 0.204 45.69 0.03 43.61 0.56 44.17 -2.08 

J0903 .4+4651 S4 0859+470 FSRQ 1.462 46.95 0.11 44.44 0.58 46.1 -2.51 

J0834.3+4221 S4 0830+42 FSRQ 0.249 45.24 0.07 43.32 0.57 44.18 -1.93 

J0920.9+4441 S4 0917+44 FSRQ 2.189 48.72 0.02 45.6 0.56 45.84 -3.12 
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J0849.0+0455 BZB J0849+0455 BLL 1.07 46.6 0.14 44.21 0.58 45.37 -2.39 

J23 19. 1-4208 PKS 2316-423 BLL 0.056 43.54 0.19 42.2 0.59 44.02 -1.34 

J2204.6+0442 PKS 2201+04 BLL 0.028 43.11 0.08 41.91 0.56 43.38 -1.2 

J2159.9+1023 CRATES J2200+1030 BLL 0.172 44.62 0.12 42.91 0.57 44.23 -1.71 


J0250.6+1713 BZQ J0250+1712 FSRQ 0.207 45.13 0.08 43.25 0.57 44.25 -1.89 

J01 12.8+3208 NRA0 62 FSRQ 0.603 46.99 0.02 44.46 0.56 45.21 -2.52 

J0909.7-0229 PKS 0907-023 FSRQ 0.957 47.13 0.04 44.55 0.57 45.3 -2.57 

J0912.9-2102 BZB J0913-2103 BLL 0.198 45.01 0.09 43.16 0.57 44.45 -1.84 

J1204.3-0711 BZB J 1204-07 10 BLL 0.185 44.91 0.08 43.1 0.57 43.84 -1.81 
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Table S2. Properties of the GRB sample. 
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000926 2.04 1.3 6.2 53.45 53.82 51.58 53.96 52.09 BeppoSAX 3,4,5,22 

010222 1.48 74.0 3.2 53.93 52.46 49.65 54.32 50.04 BeppoSAX 3,4,5,22 

011211 2.14 51.0 6.4 52.83 51.62 49.41 53.08 49.66 BeppoSAX 3,4,5,22 

020405 0.7 40.0 7.8 52.86 51.49 49.45 53.63 50.22 BeppoSAX 3,4,5,22 

020813 1.25 89.0 3.1 53.89 52.29 49.47 54.31 49.89 HETE 3,4,5,22 


021004 2.32 52.4 12.7 52.75 51.55 49.93 53.89 51.07 HETE 3,4,5,22 

031203 0.1 40.0 9.0 50.99 49.43 47.52 51.14 47.67 Integral 9,6 

030329 0.17 22.3 5.1 52.1 50.82 48.42 52.8 49.12 HETE 7,11 

050709 0.16 0.1 17.2 50.0 51.22 49.87 51.05 50.92 BAT 12 

050820A 2.62 26.0 6.6 53.98 53.13 50.95 54.73 51.69 BAT 13,14 
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090423* 8.2 10.3 > 12 53 53 > 51.3 53.6 > 51.9 BAT 19,20 

020903 s 0.251 10 < 90 49.04 48.14 < 48.14 50.6 <49.7 HETE-2 24,25 


(a) Uncertainty of 0. 1 dex as described in the text. 

(b) Uncertainty of 0.2 dex as described in the text. 

(c) Uncertainty of 0.3 dex as derived from uncertainties in 6j and E^' 

(d) Uncertainty of 0.5 dex as described in the text. 

(e) Uncertainty of 0.54 dex as derived from uncertainties in 9j and E) 
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